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ABSTRACT 

Photometrically distinct nuclear star clusters (NSCs) are common in late-type-disk and spheroidal galaxies. 
The formation of NSCs is inevitable in the context of normal star formation in which a majority of stars form 
in clusters. A young, mass-losing cluster embedded in an isolated star-forming galaxy remains gravitationally 
bound over a period determined by its initial mass and the galactic tidal field. The cluster migrates radially 
toward the center of the galaxy and becomes integrated in the NSC if it reaches the center. The rate at which the 
NSC grows by accreting young clusters can be estimated from empirical cluster formation rates and dissolution 
times. We model cluster migration and dissolution and find that the NSCs in late-type disks and in spheroidals 
could have assembled from migrating clusters. The resulting stellar nucleus contains a small fraction of the 
stellar mass of the galaxy; this fraction is sensitive to the high-mass truncation of the initial cluster mass 
function (ICMF). The resulting NSC masses are consistent with the observed values, but generically, the final 
NSCs are surrounded by a spatially more extended excess over the inward-extrapolated exponential (or Sersic) 
law of the outer galaxy. We suggest that the excess can be related to the pseudobulge phenomenon in disks, 
though perhaps not all of the pseudobulge mass assembles this way. Comparison with observed NSC masses 
can be used to constrain the truncation mass scale of the ICMF and the fraction of clusters suffering prompt 
dissolution. We infer truncation mass scales of < 1O 6 M (> 10 5 M©) without (with 90%) prompt dissolution. 
Since the NSC assembly is collisionless and non-dissipative, no relation to the process responsible for central 
black hole assembly in more massive galaxies is expected. 

Subject headings: galaxies: evolution — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
star clusters — galaxies: structure 



1. INTRODUCTION 

Photometrically distinct nuclear star clusters (NSCs) are 
ubiquitous in dynamically primitive galaxies, which are the 
galaxies that lack what can be interpreted as the structural sig- 
natures of major — "wet" or "dry" — mergers: surface bright- 
ness profiles with de Vaucouleurs-like or higher Sersic]( |T963) > 
indices, and, absent core scouring by binary massive black 
holes, 



po wer-law central luminosity cusps (e.g. 
2009a|b 



Hopkins et| 
Imaging with 



al. 2008 2009a|b and references therein), 
the Hubble Space Telescope (HST) has revealed that 75% 
of late-type (Sc-Sd) disk galaxies contain compact, lumi- 
nous star clusters at their centers with masses in the range 
8 x 10 5 - 6 x 1O 7 M , effective radii in the range l-9pc 
with a median of ~ 3.5 p c, and velocity dispersions in the 
range 13-3 4kms"' (e.g., [Boker et al.||2002[ 



in galaxies that have not experienced major mergers in the 
epoch follo wing the initial burst of star formation (see, e.g., 
|Koda et aL||2009| l, but have a well-defined dynamical cen- 
ter. The faintest dwarf spheroidal satellites of the Milk y Way, 
with th e exception of the Sagittarius Dwarf (e.g., |Bellazzini et 
|al.|2008[ and references therein), do not contain NSCs, even 
though they do seem to form a continuous structural family 
of stellar systems with the more massi ve dwarf sphe roidals 
that d o contain NSCs (e.g., Kor mendy et al.||2009| [Wyse 



2010 and references therein). While the spheroidal galaxies 
in Virgo are typically non-rotating (Geha et al. 2003), their 



2004 Walcher 



|et al.||2005] ). The clusters 5 ages range widely from 4x10' 
to 10 lu yr, the metallicities average at (Z) = 0.015 with a sig- 
nificant scatter, the average star formation rate over the last 
100 Myr is (S SFR ) = 2 x 10" 3 M Q yr _1 , and older stars dom- 
inate the cluster mass. Similarly, an HST survey of stellar 
nuclei in spheroidal galaxies in the Virgo and Fornax clusters 
( |Ferrarese efaL] 2006a; C ote et al.|2006||Miller & Lotz|2007"] l 
has revealed stellar nuclei th at on a verage cont ain ~ 0.2% 
of the ga laxy mass ( Ferrarese et al.|[2006b Wehner & Har- 



nucleus-subtracted surface brightness profiles, with Sersic in 
dices n ~ 1-2, resemble those of disks. This motivates a 
unified approach to addressing NSC formation in these two 
morphological classes. In disks, an excess surface bright- 
ness above the inward-extrapolated exponential law of the 
outer disk occurring at radii ~ 100-500p c, variously char- 
acterized as "c entral light excess" (|Boker et al.||2003|l or a 
"pseudobulge" (Kormendy & Kennicutt 2004; Fisher & Drory 



ns 



2006) a nd have a m edian effective radius of 4.2 pc (see 
aIsb, |Geha et dl2002) . Also, HST imaging of nearby low- 



luminosity dwarf galaxies with absolute magnitudes -16 < 
My < -13 ( [Georgiev et al.||2009| l has revealed NSCs with 
masses ~ (10 b - 1O 7 ).M0, with a majority clustering at the 
lower end of the mass range, and half-light radii 1 .5 — 4.7pc. 
It seems that photometrically distinct NSCs are common 



[20081 [20101 IFisher et al.|2009||Weinzirl et al |2009| l, may pro- 
vide additional clues, because the process that produced the 
NSC could also contribute to the formation of the pseudob- 
ulge. Pseudobulges can be distinguished from classical bulges 
through the ir cold, rapidly rotating k inematics and low Sersic 
index (see [K ormendy & Fisher 2008 and references therein). 
Kormendy et al. (2010) find that the pseudobulges in the sam- 
ple of pure-disk (Scd) galaxies contain < 3% of the galactic 
stellar mass, which is smaller than the stellar mass fraction in 
classical bulges. 

Here we discuss the formation of NSCs in the context of 
normal star formation in which most stars form in clusters. 
A successful theory of NSC formation should explain the ob- 
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served uniformity of their properties. The smallness of their 
masses in comparison with the total stellar masses of the host 
galaxies is a "small parameter" characterizing the assembly of 
cosmic structure that requires explanation. A theory of NSC 
formation could help constrain the earliest stages of the host 
galaxy's formation. NSCs are interesting because they seem 
to betray either an inefficiency in the transport of baryons 
(gas, stars) to the center of the galaxy, or perhaps an ineffi- 
ciency of star formation once gas has arrived at the center of 
the galaxy. The NSC phenomenon may be related to the rel- 
atively slow, or "secular" buildup of the host galaxy's stellar 
mass, since the galaxies that form the bul k of their stars on dy- 
namica l time scales, as in mergers (e.g., Hopkins et al .|2008 



raling star clusters. In a similar spirit but different context and 



2009a ) and cold accretion from the intergalactic environment 



(e.g., |Dekel et al.|[2009al |Ceverino et~aLlj20T0] >, should turn 
out to be more concentrated than the NSC host galaxies and 
should contain "classical" bulges structurally equivalent with 
ellipticals. 

NSCs can grow from the local interstellar medium if the 
galactic gas acc retes and accumulates at the center (e.g., 
Seth et aL]|2006| ). In gas disks, sufficient gas accretion to 
form NSCs may simply be driven by magnetic stresses in 
the gas that are ampl ified by the magnetorotational instability 
(Milosavljevic 2004), w hich is generic in differentially rotat- 



ing galactic disks (e.g., |Kim et aT]|2003 ). If the inner bary 



onic disk is globally self-gravitating, accretion can be driven 
by torques associated with stellar and gaseou s bars (e.g., Ko- 
|rmendy & Kennicutt||2004| |Schinnerer et al. 2007] and ref- 
erences therein). NSCs in late-type disks typically contain 
young stellar comp onents, and their star formation se ems to 
be intermittent ( |Rossa et al.|2006[|Walcher et al.|2006| ), which 
is consistent with the gas accretion scenario. Cosmological 
hydrodynamical simulations of dwarf and late-type galaxy 
formation com monly yield nucleus- (or "bulge"-)disk struc- 
tur es (e.g.,|Govemato et aIlpO^|2Wl|2009j|Kaufmann et| 
|aIJ|2007T|Brooks et al.||2009| |Ceverino et al.||2~010| ), where 
the nucleus forms from the gas that self-gravitational torquing 
has driven into the center of the galaxy, whereas the disk 
forms from the gas that manages to stably circularize in the 
combined gravitational potential of the dark matter halo and 
gaseous stellar nucleus. The resulting nuclei normally con- 
tain larger fractions of galaxy mass than the observed nuclei 
in late-type disks and spheroidals. Frequently, a gas expulsion 
by feedback from star formati on is invoked to explain t his ap- 
parent discrepancy (e.g., Zavala et al.|20 08 ; Dutt on|2009[ ). 

Here we do not study gas accretion, but instead investigate 
the possibility, which is necessitated by the clustered charac- 
ter of star formation, that a substantial fraction of the NSC 
mass could have assembled in off-center stellar associations 
(clusters) that subsequently migrated into the galactic cen- 



ters and merged into an NSC as intact en tities (e.g., Tremaine 
et al.||1975||Lotz et al.||2001| |Capuzzo-Dolcetta & MioccrT 



2008ajbj). The assembly of NSCs through the merging of 
clusters is consistent with the observation that nuclear cluster 
phase space densities are on average somew hat smaller than 
those of globular clusters ( Walcher et al. 2005 ) because phase 
space density tends to increase in collisionless mergers. The 
NSCs thus assembled should inherit the orbital angular mo- 
mentum of the migrating disk clusters and should be rotating 
(s ee, e.g |Seth et al.|2008| . 

Belda] ( |2010| ) carried out A^-body simulations of star clus- 
ter orbital decay in the background of field stars in a disk 
galaxy embedded in a dark matter halo. His results suggest 
that NSCs could have formed from stars delivered by inspi- 



follow i ng an earlier p roposal by Noguchi ( 1999 ), Immeli et al. 
d2004j), [Bou rnaud et al. (2007), Elmegr een et al'| ( |2008| ), and" 
Cever ino et al.| (2010) carried out numerical simulations to 
find that the migration and central merging of giant star clus- 
ters and gravitationally bound gas clumps in massive, rapidly- 
star-forming disks can give rise to classical-bulge-like, dy- 
namically hot central stellar systems. If giant cluster migra- 
tion can produce classical bulges in gravitationally unstable 
disks that are fed by a relatively smooth accretion from the 
extragalacti c environme nt — and not by major mergers (see, 
e.g., |Genzel et al.|20" 08 ) — then this warrants an investigation 
of cluster migration in lower-surface-density galaxies, in hope 
that an improved understanding of their morphological trans- 
formation can be gained to elucidate why these galaxies end 
up lacking classical bulges and remain rotationally-dominated 
and dynamically cold. 

To this end, we develop a crude toy model to calculate the 
influx of migrating young clusters formed in a galaxy that has 
experienced an instantaneous star formation episode, but our 
results are relevant for galaxies, like the prototypical late-type 
disk galaxy M33, in which star formation has been steady and 
ongoing over the life of the galaxy, although the timescale 
for the accumulation of NSC mass will differ from the one 
estimated here because of the dependence, in this case, on 
the star formation rate in the disk. We assume that all stars 
form in clusters and that the clusters eventually get disrupted 
in the galactic tidal field. We further dis cuss the assumptions 
and limitations of our model in Section |3.4| Gas expulsion, 
stellar mass loss, dynamical evaporation, and external pertur- 
bations lead to eith er rapid or progressi ve tidal mass loss in 
the clusters (e.g., [Fall & Z hang 2001 ; Baumgardt & Makino 
[2003| pVIcMillan & Portegies Zwart|2003[|Gieles et aTp 006a). 
If a cluster is massive, it can migrate radially inward prior 
to complete disruption. For example, we will find that for 
our galactic models and parameters (see Sections [3] and |4| 
a 10 4 M Q cluster never migrates more than lOpc; a 10 5 M Q 
cluster can migrate and reach the center from ~ few x lOpc, 
a 10 6 M Q cluster can do so from ~ 100pc and a 10 7 M© clus- 
ter can reach the center from > lkpc. The cluster's stars are 
then deposited at galactocentric radii traversed by the migrat- 
ing cluster. Since even in very late-type disks (Hubble type 
Scd) circular test particle orbits rotate differentially through- 
out, cluster disruption is inevitable except if the cluster forms 
at, or migrates into, the galactic center, where tidal shear pre- 
serves the symmetries of an axisymmetric cluster. Clusters 
with off-center birthplaces can avoid complete disruption if 
they migrate to the dynamical center quickly enough. The 
progressive tidal disruption of migrating clusters before they 
have reached the galactic center transports mass radially in- 
ward; we propose that this process is a driver of pseudobulge 
growth in disks. 

This work is organized as follows. In Sectionj2]we provide 
an empirically and theoretically motivated description of clus- 
ter formation, migration, and dissolution. In Section[3]we in- 
troduce our models for a spheroidal galaxy and a disk galaxy 
initially lacking an NSC, and estimate the effects of galactic 
stellar mass redistribution due to cluster migration. In Section 
Hlwe show our results for representative galaxy models. In 
Section [5] we compare our results with observations of NSCs 
and pseudobulges and provide a brief discussion of implica- 
tions for the pseudobulge and classical bulge dichotomy and 
for massive black holes, and in Section [6] we summarize our 
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main conclusions. 

2. CLUSTER FORMATION, MIGRATION AND DISSOLUTION 

We assume that all stars in our model galaxies formed in 
clusters, and introduce a framework based on previous em- 
pirical and theoretical work for treating the formation, migra- 
tion and dissolution of these young stellar clusters. In Section 
|2.2| we discuss a theoretically motivated description of clus- 
ter migration in spheroidal and disk galaxies, taking into ac- 
count the dominant form of dynamical friction torque ac ting 
on a migrating cluster in each type of galaxy. In Section 2.3 
we describe an empirically motivated model of cluster disrup- 
tion based on o bservations of young clusters in disk galaxies. 
In Section |2~T| we introduce the initial cluster mass function, 
which seems to be generic to star formation in all galaxies, 
including the Milky Way. 

2.1. The Initial Cluster Mass Function 

The initial cluster mass function (ICMF) contains about 
equal mass on all cluster mass scales, i.e., 



the column density N ~ 10 cm that is required for the pres- 



dn 
~dM 



oc M a for M min <M <M n 



(1) 



with a ~ 2 (e.g 
[eTaLl[2003t |Rafe 



Bik et al. 2003 
ski & Zarits" 



03l|deGrijset al.|2003| 
kfpnU5n Chandar STST 



Hunter 
2010), 



where M denotes initial cluster mass and M m i n and M„ 
galaxy-dependent cutoffs. We adopt a = 2 and M m \ n = 100 M Q 
in what follows; our results will depend only weakly on M m j n 
since clusters with masses below ~ 10 4 M Q are disrupted be- 
fore they can migrate to the center from larger radii and, if 
they do reach the center, contribute an insignificant fraction 
of the total NSC mass. The definition of a high-mass cutoff of 
the ICMF and the statistical significance of observational evi- 
dence for such a cutoff have been the subject of debate. We as- 
sume here that the ICMF is indeed subject to high mass trun- 
cation, and treat M max as a free parameter that varies over the 
range 1O 4 M - 10 7 M Q , which includes the cutoffs reported 
for the nearby galaxies. We will find in Section|4]that masses 
of nuclear clusters in our calculations are most sensitive to 
the truncation mass, and thus the photometry, of spheroidal 
and disk galaxies. If combined with a theoretical model of 
cluster migration and dissolution, this sensitivity can be uti- 
lized to indirectly constrain the ICMF truncation mass scale 
in these galaxies. 

The maximum cluster mass forming in NGC 6946, M51, 
and the Antenna e has been estimated f rom observations to be 
10 6 M o (|Gieles et al.|2006b|cl. Frequently, an ICMF 



with an exponential cutoff dn/dM oc M~ a e~ M l M * is found to 
adequately approximate the truncation of the ICMF at the 
high-mass end. In spirals and irregulars the mass truncation 
scales M* and M max are in the range ~2x 10 s - 10 7 ^tfp 
(see, e.g., [Z hang & Fall|1999||Bik et al.|2003)|de Grijs et al. 
12003) |Dowell et al.|2008[|Larsen|2009[|Gieles|2009| >. The ex- 
ponential mass truncation scale, which has typical values of 
M+ ~2x 1O 5 M , is larger in denser starburst environments 
( |Larsen|2009| l. In the present work, we adopt the abrupt trun- 
cation with M max , but we do not anticipate that the results 
would be substantially different if exponential truncation were 
used. 

The maximum mass of clusters that formed in the assembly 
of nucleated spheroidals is unknown, but a brief theoretical 
speculation may be in order. In a gas disk with surface density 
E gas close to the critical value < 



ence of a self-shielding cold phase (see, e.g., Schaye||2004 



and Toomre parameters Q < 1, the truncation scale M u 
M*) should scale with the Jeans mass of the disk, 



(or 



Mrr 



r /47r 2 G£ gas 

VSF 



^4xl0 6 M o / SFr 



VlOMopc" 



K 

io- 15 f 



(.2) 



where /sf = 0.1/sf-i is the star forma tion efficiency and k 
is the epicyclic frequency (see, e.g., Elmegr een et al.[ 2008). 
The reference epicyclic frequency in equation was se- 
lected to correspond to the characteristic average total mass 
density ~ 0.1M Q pc~ 3 in the inner 300pc of the (largely non- 



nucleated) dwarf spheroidal satellites of the Milky Way ( Stri 
gari et aL|2008) , but of course, the proto-spheroidals would 
have been characterized by larger, radius-dependent values 
of k. If the gas disk mass is a fixed proportion of the to- 
tal (e.g., dark matter) enclosed mass, and the vertical scale 

height of the disk is h, then k oc E ga s /T 1 ' 2 , which would im- 
ply a weaker dependence M max oc £ gas h 2 on the gas surface 
density. It is plausible that proto-spheroidals assembled from 
an ICMF reaching the cluster mass scale associated with the 
gas disk Jeans mass estimated in equation (|2]). 

2.2. Cluster Migration 

A gravitationally-bound star cluster migrates on a time 
scale 

dJ 



d\nR 



IDF 



(3) 



where R is a characteristic size of the orbit — the equivalent 
of the semimajor axis in a non-Keplerian potential, J is the 
angular momentum of the cluster, and (Tdf) is the orbit- 
averaged dynamical friction torque. All galactic mass compo- 
nents (dark matter, gas disk, stars) respond dynamically to the 
cluster but the mechanism of torque coupling varies. We sep- 
arately consider the dynamical friction torque in spheroidal 
and disky mass components. 

In spheroidal stellar systems or dark matter halos, A^-body 
simulations have shown that the torque from a nonrotat- 
ing collisionless halo can be heuristically described with the 
Chandrasekhar formula for dynamical friction, 



'DF,Chandra : 



4tt ln(A)x( V)G 2 M 2 p 



rxV, 



(4) 



10M Q pc 1 corresponding to 



where p is the local combined density of stars and dark matter, 
|V| ~ flR is the velocity of the cluster, ln(A) is the Coulomb 
logarithm, and x00 is the mass fraction of stars or dark matter 
particles with velocities less than V. Because the kinematic 
structure of actual halos differs from the premises of Chan- 
drasekhar's derivation, A^-body simulations are necessary to 
obtain the correct normalization of the torque amplitude. The 
numerically evaluated torques ( e.g., [Velazquez & White|1999{ 
IPenarrubia et~ai~||2lX)2"l [2004; Spi nnato et al.||2003| > caiTbe 
modeled with Chandrasekhar's formula if the Coulomb log- 
arithm is treated as an empirically-calibrated free parameter. 
The values of the Coulomb logarithm obtained from these Af- 
body simulations are ln(A) ~ 2— 7. 

In differentially rotating gas or stellar disks, the torque is 
provided by the flow in the corotation region and by angu- 
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lar momentum transfer at Lindblad resonances (e .g., Goldre 
|ich & Tremaine|1980|jQuinn & Goodman|1 986). Disk clus 
ters with masses M < 10° M Q have Roche tidal radii r t ~ 
[GM/(dil/dlnR) 2 ]^ 3 that are smaller than the thickness of 
the disk. In this regime, the disk torque is a generalization of 
the "Type I" torque acting on small planets in protoplanetary 
disks (e.g., Tanaka et al. 2002; D' Angelo et al.|200 3 1 Baruteau 
|& Mas set 2008 ) to non-Keplerian disks and is proportional to 



TbF.disk OC 



G 2 M 2 S 



(5) 



where £ is the surface density of the gas or stellar disk, and 
a is the gas sound speed or the stellar velocity dispersion in 
the disk. In what follows we take the constant of proportion- 
ality in equation |5]l to be equal to unity. This choice seems 
consistent with the numerical cal ibration of the torque in a 
collisionless particle disk by Donner & Sundelius ( 1993 1 and 



Wahde et al. ( 1996 1 if the gravitational softening employed in 
their calculations is equated with the vertical thickness of the 
disk. 



Finally, we note that |Milosavljevic| ([2004 ) argued against 
the migratory scenario in disks on the grounds of long migra- 
tion times. There, migration from a distant location in the disk 
was envisioned, and the treatment ignored the enhancement in 
dynamical friction force due to the stellar accumulation of ear- 
lier migrating clusters in the central kiloparsec of the galaxy. 
In the present picture, clusters migrate from a range of radii 
in the galactic disk, but only those that form closest to the 
galactic center reach the center and merge with the nuclear 
cluster. 

2.3. Cluster Dissolution 

The timescale on which a cluster is disrupted in the galac- 
tic tidal field has been determined empirically by modeling 
the luminosity and age functions of clus ters in nearby disk 
galaxies (see, e.g., La mers et al.|2005a|bj ). Theoretical mod- 
els tracking stellar and dynamical evol ution of a cluster in a 
tidal field (e.g., Gieles et al.|2008"l|Gieles & BaumgardtpOW 



have reproduced the observationally inferred variation of the 
dissolution time with cluster mass, 



f dis - to 



M 



(6) 



where 7 ~ 0.62 (Boutloukos & Lamers 2003 de Grijs & An 



ders[2006). The normalization to varies between galaxies, 
which seems to be a consequence of the variation in the tidal 
field strength and of cluste r-scale density inho mogeneities of 
the galactic environment. |Gieles et al. ( |2008| l assessed the 
role of the tidal field by comparing cluster lifetimes in sev- 
eral galaxies with the inverse angular frequency of circular 
orbits at observed cluster radii and found consistency with 
the linear relation t^ oc fT 1 . The residual variation of fdis 
between galaxies has been as cribed to disruption by giant 
molecular clouds ( |Gieles et al.|2006al and references therein) 
which we ignore. The tidal radius really depends on the de- 
gree of differential rotation , r t = [GM /{dVL/dXnR) 2 } 1 ^ (e.g., 
IQuinn & Goodmanl|1986|l, and thus one would expect that 
tp oc \dfl/dlnR\~ l . |Lamers et al. ( 2005a 1 modeled the ob- 
served cluster population in M33 and found, at an arbitrary 
radius, t « f&s/tt(R) with / dis ,M33 ~ 0.16. 

The coefficient f& s , which can be separately inferred from 
the observed young cluster populations and from theoretical 



calculations, encapsulates the detailed mass loss mechanics 
in a cluster embedded in a galactic tidal field and evolving 
through internal two-body relaxation. From the theoretical 
viewpoint, the coefficient depends on the initial stellar mass 
function, on the initial cluster structure, and on the (possi- 
bly time-dependent) gravitational potential within and near 
the cluster. Under idealized assumptions, the coefficient has 
been estimated with jV-body s imulations (Portegies Zwart et 
|al.l 1998| [2002} [Baumgardt & Makino 2003). The theoretical 
estimate /dis ~ 0^ in|Lamers et al.|(|20 05a), which is based 
on the results of Baumgardt & Makino (2003 ), is larger than 
the empirical value. The clusters originating in M33's cen- 
tral disk could be denser an d more resistant to tidal disruption 
than those in the sample of Lamers et al. ( 2005a I. 

The above relations were obtained for clusters originating 
in a galactic disk. For lack of an equivalent empirical result 
for spheroidal starbursts, we assume that these relations hold 
universally and adopt the crude relation 



to = fdi 



dQ 



d\nr 



(7) 



where we take f^ s = 0.2. We have compared our simpli- 
fied form of the dissolution timescale defined in equations 
^ and Q for a cluster with initial mass M = 10 s M Q and 
angular velocity fi(r) = [G ; M ga i i o(r)/r 3 ] 1 ' /2 , where M g!k \fi(r) = 
M+fl(r) +Mdm(>") is the total initial galactic mass within a 
sphere of radius r (see Section 131, with the timescale from 

_ n 1 1 . 1 M . a _ . . .. rr 



Lamers et al.| d2005a|bj>, where t' dis = t' 4 (M/10 4 M Q ^ (Bout- 



loukos & Lamers|2003| l, t' A = 1.355 yrx lO 4 ^" 1 ^/! yrp* 77 
( |Lamers et al.|2005b|l , and fg = /di s f2(r) -1 , and found it to be 
consistent to within 20% in the range of radii where we find 
that significant cluster migration takes place. 

3. INITIAL AND REDISTRIBUTED GALACTIC STELLAR MASS 

Here we describe our models for the initial density profile 
of a sph eroidal galaxy (Section [3. \\ and a disk galaxy (Sec- 
tion 3.2 1. The initial density lacks a stellar nucleus; its inner- 
most baryonic density profi le is an extrapolation from larger 
radii. Then in Section 13.31 we describe our method of dis- 



tributing initial clusters consistent with these surface density 
profiles and modeling cluster mass redistribution during mi- 
gration and dissolution. We include the growing central clus- 
ter mass in the mass distribution affecting subsequent migrat- 
ing clusters, but we ignore the response of the dark matter 
halo to the baryonic collapse and subsequent mass redistribu- 
tion by cluster migration. 

3.1. Spheroidal Galaxy Model 
Photometry of spheroidal galaxies in the nearby universe 



(Cote et al. |2006 Ferrarese et al.||2~006a| [Kormendy et al. 
2009 ), combined with the fact that in these relatively old stel- 



009 ), combined with the fact that in these relatively old ste 
lar systems the stellar mass-to-light ratios do not vary signifi- 
cantly within galaxies, has shown that their mass surface den- 
sity profiles are well desc ribed by a two-co mponent model 
( Balc ells et al.|2003| |Graham & Guzman|2003l 



£*W = S S ersic(^) + S nuc l(/?). 



(8) 



The Sersic law component T,s ers i c (R) ^ ex P[ - 0^/^s) 1// "] 
HCaon et al.TT993l ,where n is the Sersic index and R s is a scale 
radius, contains most of the mass, and the com pact nuclear 
component £ nuc i(j?) contains a small fraction ( Ferrares e et al.| 
2006b; [Wehner & Harris||2006[ ) of the mass. In spheroidals 
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the Sersic index varies in the range n ~ 1 - 2 ( Fe rrarese et al. 
2006aj Kormendy et al.|2009 1. We work within the paradigm 



in which a stellar nucleus is not present prior to the forma- 
tion of most of the galactic stars. Thus, we consider an initial 
galaxy with a spherically-averaged stellar density profile that 
lacks a nucleus. We employ an approximation to the depro- 



jected Sersic profile i n the form (Prugniel & Simien 1997 
|LimaNetoetal.|1999) 



exp 



i/n 



where 



p= 1-0.6097 «"'+ 0.05463 n" 



(9) 



(10) 



Here, p sp h is normalized to the total stellar mass of the galaxy 
M sph = 4tt J °° Pi , (r)r 2 dr = 47r«^p sph r[(3 - p)n], and r(a) = 
f^t^e-'dt is the gamma function. Quantities pertaining to 
the initial galaxy, before migration and disruption of clusters, 
are denoted by the subscript 0. The stellar mass within radius 
r is 

M*, (r) = 47rn^p sph7 [(3 - p)n, (r/R^ 1 "] , (11) 

where j(a, x) = ^t a ~^e~'dt is the lower incomplete Gamma 
function. 

The total initial density profile includes a dark matter com- 
ponent 

Po(r) = p+, {r) + p DM (r), (12) 

whe re we model the dark matter density profile with the NFW 
law ( |Navarro et al.|l997 i PdmM = p s /[(r/r s )(l + r/r s ) 2 ] with 
scale radius r s , virial radius r v i r , and concentration c = r wa /r s . 
Spheroidals are often dark matter dominated, especially at 
the low-luminosity end, i.e., Pdm(>") ^ P*.o(r)- Thus, for 
galaxies in which dark matter is the dominant mass compo- 
nent, the cluster migration time can be modeled with the N- 
body-simulation-calibrated Chandrasekhar dynamical friction 
torque in equation Q. 

3.2. Disk Galaxy Model 

The surface density profiles of late-type spiral disk galax- 
ies can be fit with a Sers ic-type expone ntial stellar disk and 
a nuclear component (B oker et al.|2003] ) as in equation dHJ. 
We consider "pure" disk galaxies, which do not have a stellar 
bulge or a preexisting pseudobulge; we do not consider clus- 
ter formation in the galactic halo. The initial stellar density 
profile without the nuclear component is then 



E*.o(/?) = Sdiskexp 



i/n 



(13) 



where Sdisk = A^disk/[27rn^r(2n)]. The corresponding initial 
stellar mass profile is 



M*, (/f) = 27m^£ disk r[2n, {R/R % )V n ], 



(14) 



where T(a,x) = T(a)—j(fl,x) is the upper incomplete Gamma 
function. The total initial density profile is then 



E (/f) = £*,o(*) + E DM (/f), 



(15) 



where Edm(^) is the projected NFW profile. Disk galaxies 
generally have an exponential profile with Sersic index n ~ 1 . 

The total density profiles in disk galaxies may be dominated 
either by dark matter or by luminous matter in their central 



regions. Even if dark matter dominates the spherically aver- 
aged density profile, in disk galaxies it provides a relatively 
smaller contribution to the torque if the cluster is embedded 
in the disk. This is because the torque coupling of the migrat- 
ing star cluster with the flattened disk, equation i s more 
significant than that with the dark matter halo (see, e.g., |Bekki] 
20 10) . We use the disk torque to model the migration time in 
disk galaxies. 

3.3. Stellar Mass Redistribution by Cluster Migration 

We assume that all stars form i n cl usters with masses be- 
tween M m ; n and M max (see Section |2Tlj i. Let d 2 n /dMdr be the 
number of clusters per unit initial cluster mass per unit radius 
that have formed in a galaxy such thatf °° (d 2 n/dMdr) dr is 
proportional to the ICMF in equation ffl, while 



M m 



A/,„ 



d 2 n 
dMdr' 



MdM- 



(4Trr 2 p+x)(r), spheroid, 
\ 27rr£*,o(r), disk. 



(16) 



We sample cluster masses and initial radii randomly ac- 
cording to the distribution d 2 n/dMdr. The total num- 
ber of clusters in our simulations is J M m " (dn j dM) dM w 
M tot /[M min ln(M max /M min )], where M tot is the total stellar 
mass, M sp h or Mdi s k, of the galaxy. At first, we do not ex- 
plicitly take into account the possible prompt, possibly mass- 
independent (for a discussion of mass-dependent prompt dis- 
solution see, e.g., [Goodwin 2009) dissolution of clusters ("in- 
fant mortality," see, e.g.,|Fall et al.|2 005 ; Bastian & Goodwin 
2006)|Chandar et al.|2006HGoodwin & Bastian|2006 t IGieles 
et al.|2007MPellerin et al.|2007HWhitmore et al.|2007 
gardt et al.|2008||de Grijs & Goodwin|2 0081 |Lamers 



2008; I 
•ompt di 



Baum- 



2009). 



Then, we briefly assess the effects of prompt dissolution . 

We proceed to model cluster migration and dissolution by 
following the cluster orbital decay 



and mass loss 



dr 
dt ~ 

dM 
dt 



'mig 

M 

fdis ' 



(17) 



(18) 



where f m j g and t<a s are, respectively, the migration time scale 
in equation ([3]l and dissolution time scale in equation |6|). 

We calculate the orbital decay of clusters in the order of 
increasing migration time (see equation |3)) evaluated before 
migration has occurred. We keep track of the mass tidally 
stripped from the cluster with kth shortest migration time by 
calculating its contribution to the stellar density profile 



Anr 2 



dM k 
dt 



1 dr k 
dt 



(rk < r k ,o), 



(19) 



where r^o is the cluster formation radius and k = 1,2,.. .n tot . 
We assume that a cluster has been fully disrupted when its 
mass falls below M m ; n = 100 M Q since in clusters of this size 
?mig ^ fdis and thus the cluster is not able to migrate sub- 
stantially, if at all, before disruption; this residual mass is 
deposited at the radius of disruption. Some migrating clus- 
ters avoid complete disruption until they reach the innermost 
radius of our grid, r m ; n = 1 pc. The migrating cluster mass 
that reaches the central parsec is added to the central point 
mass M cent . To account for the influence of the growing nu- 
clear component, including M cent and the mass deposited at 
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larger radii, e.g., at r < 100 pc, on the migration and disrup- 
tion time scales of more slowly inspiraling clusters, we cal- 
culate the tidally stripped mass of faster migrating star clus- 
ters with smaller f m j g first. After the migration and disrup- 
tion of k clusters has been computed, we sum up their contri- 
butions to the modified stellar density profile of the galaxy 
Ap^ = Ap^fc/, and set the stellar density profile in 

which the (k+ l)st cluster migrates to max(p+ q, so that 

the late-migrating clusters do so in a galaxy modified by the 
central concentration increase from the early-migrating clus- 
ters. We take the final density profile of the galaxy, after all 
clusters have been disrupted, to be p+ = Y^k^P*,k- This pre- 
scription does not account for the escape of the material liber- 
ated by stellar mass loss from the galaxy. To model the effect 
of this escape, we repeat the calculation after requiring that 
the first 50% of the mass to be stripped from a cluster leave 
the galaxy. 

We treat the mass delivered to the central parsec as an "un- 
resolved" stellar nucleus component with mass M cen t. We 
model its surface density profile with the function 



S C ent(^) — 



(20) 



where R cenl is the half-light radius of the unresolved compo- 
nent. We set this radius to R cent = 2 pc, which is chosen to 
be larger than our innermost grid radius but smaller than typi- 
cal observed hal f-light radii of NSCs, e.g., R nuc > lOpc (e.g., 
|Geha et al.|2002| >. We add this unresolved component to the 
surface density profile calculated from the mass deposited by 
migrating clusters at r > 1 pc, so that 



S u „,ai(R)= / pA(R 2 +z 2 ) 1/2 ]dz+Z ce JR). 



(21) 



Our surface density profiles do not take into account any 
smearing by the point spread function (PSF). If the PSF width 
is comparable to the resulting NSC radius, then smearing by 
the PSF must be taken into account when comparing synthe- 
sized surface brightness profiles with the observed ones. Be- 
cause our results are affected by a variety of crude approxi- 
mations we do not proceed to model the effect of the PSF. 

3.4. Limitations of the Model 

Before proceeding to discuss our results we would like to 
highlight the assumptions and limitations of our model. We 
have assumed that all clusters form simultaneously and thus 
we calculated the orbital decay of the clusters in order from 
fastest-migrating (as calculated before any migration has oc- 
curred) to slowest-migrating. This assumption is applicable to 
those dwarf galaxies that consume gas and form stars in a sin- 
gle star formation episode. In other galaxies, including late- 
type disks, however, cluster formation is ongoing (e.g., in the 
galaxy M33), or is triggered by galactic mergers (e.g., in the 
Antennae galaxies). Our model could be adapted to approxi- 
mate ongoing star formation by calculating cluster orbital de- 
cay and the associated mass transport in random order so that 
an older, less massive, and slower-migrating cluster can mi- 
grate and deposit its mass before a younger, more massive, 
and faster- mi grating cluster does. Without migration-time or- 
dering, our resulting NSC masses are similar to those in the 
corresponding ordered models. We do not attempt to estimate 
the time scale for the build-up of NSC mass in this scenario 
because the total timescale, although certainly longer than in 



the instantaneous star formation scenario, would depend on 
additional parameters — the star formation rate and the rate of 
gas accretion onto the galaxy — that are bound to vary between 
galaxies. 

In an effort to focus on the variation in the galactic surface 
density profile and NSC mass buildup as a function of the 
ICMF truncation mass scale M max , we do not study their de- 
pendence on other parameters such as galaxy mass and initial 
concentration (Sersic index), or on parameters, such as the 
coefficients characterizing the amplitudes of the dynamical 
friction torques driving cluster migration, that are subject to 
theoretical uncertainty. In the same spirit, we do not consider 
star formation episodes resulting from galactic mergers as this 
would require separate tracking of the galactic stellar and gas 
masses. The role of gaseous accretion and galactic mergers in 
galaxy assembly is best addressed with co mprehensive cos 
mological hydrodynamical simulations (e .g., Governat o et al. 
20041 [20071 [2009] [Kaufmann et al. 2007; Brooks et al. 200' 



Ceverino et al.|2010 1, as analytical and semi-analytical mod- 
els are inadequate in this context. We also do not model the 
response of the dark matter halo to the initial baryonic infall 
and to the subsequent evolution due to young cluster migra- 
tion and stellar mass loss. Our model for progressive clus- 
ter migration and decay is crude and thus we will compare 
only the general characteristics of the resulting surface den- 
sity profiles with those of observed galactic profiles. We cau- 
tion against comparison with the detailed features of observed 
galaxy photometry. 

4. RESULTS 

We are ready to define the parameters of our sample cal- 
culations and provide results for model spheroidal and disk 
galaxies forming at high and low redshifts. For a fixed ef- 
fective radius of the initial stellar density profile, the forma- 
tion redshift determines the average density of dark matter 
within the effective radius. High-redshift models have dark 
matter halos with lower concentration and smaller virial ra- 
dius than low-redshift models. The parameters Mh a i , c, and 
r v i r , in general, depend on the assumed cosmological model 
and on the redshift of halo formation. We assume a flat uni- 
verse : with Q m = 0.274, Q A = 0.726, and h = 0.705 ( [Komatsu| 
|et al.|2009) . The mass of the halo is related to the virial radius 
by Mhaio = 57i"A c (z)p crit ^ ir , where p a -it and A c (z) are, respec- 
tively, the critical density of the universe and the mean over- 



density within the virial radius at redshift z (Bryan & Nor- 
man 1998). For example, A c (0) « 100 and A c (5-10) w 177. 
The median concentration parameter for all high-redshift dark 
matter halos is c ~ 3.5, but at low redshift varies steeply with 
halo mass (e.g., |Bullock et al.||200l[ |Zhao et al.||2003| [Gao] 
|et al.||2008[ [Klypin et al.||2010[ 
models are for z = 6 and z 



Our high- and low-redshift 
Irrespectively. 
We pr esen t o ur re sults for spheroidal and disk gala xies in 
Sections 4.1 and 4.2 respectively. Then in Section 4.3 we ex- 
plore the observationally-motivated scenario in which a large 
fraction of the clusters dissolve promptly following forma- 
tion. 

4.1. Spheroidal Galaxies 

In Figure [TJ left and middle panel, we show the projected 
surface density profiles for stars before migration and disrup- 
tion (M sph = 1O 9 M , R s = 0.5 kpc, n = 1.5), for dark matter 
(M halo = 1O 1O M , high-z: c = 3.5, r vir = 10 kpc; low-z: c= 15, 
r v ir = 60 kpc), and for stars after migration and disruption, as- 
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FIG. 1 . — Projected initial surface density profiles of stars (thick solid line) and dark matter (thin solid line) for a proto-spheroidal galaxy and the final stellar 
surface density profile of the galaxy after all the migrating clusters have been disrupte d, for ICMF high-mass truncation of M max = 10 4 Mq (dashed line) , lQ 5 Mp, 
(dash-dotted line), 10 6 Mq (dotted line), 10 7 Mq (dash-triple-dotted line; see SectionlzTJ, and the effective Coulomb logarithm InA = 5 (see Section |2!2) . We 
show spheroidal galaxies formed at hig h red shift (left panel) and at low redshift (middle panel) ignoring galactic mass loss due to stellar evolution, and assuming 
50% mass loss (right panel; see Section[53J. 



suming In A = 5 and M max = (10 4 , 10 5 , 10 6 , 1O 7 )M . To sim- 
plify our calculation, we take p = 2/3 instead of the value 
p ss 0.62 implied by equation (JTOj. In these models we have 
not taken into account galactic mass loss resulting from stellar 
evolution. To explore the consequences of potential galactic 
loss of gas liberated by supernovae and AGB stars, in Figure 
[T[ right panel, we repeat the calculation for the high-redshift 
model assuming that the first 50% of the mass dissociated 
from each migrating cluster is lost, i.e., leaves the galaxy. The 
models shown in Figure[T]do not allow for prompt cluster dis- 
solution; we explore the impact of prompt dissolution in Sec- 
tion [43jbelow. 

The final structure of the galaxy, after all clusters have mi- 
grated and dissolved, is a strong function of the ICMF trunca- 
tion mass scale M max . For M max = 10 4 M Q , the excess surface 
density above the initial stellar density profile is present only 
in the central ~ lOpc and is small. At the other end of the 
range of ICMF truncation mass scales, forM max = 1O 7 M , the 
surface density profile of the entire galaxy has become slightly 
more concentrated, but still approximately follows the initial 
Sersic profile at large radii. ForM max = (10 4 - 10 7 )M max , the 
resulting surface density excess after cluster migration and 
disruption appears as a distinct departure from the inward- 
extrapolated outer Sersic law at R ~ 20-200pc in the high- 
Z model and R ~ 30-300pc in the low-z model; the ex- 
cess surface density tota i - S cent approximates power laws 
cx R~°- 6 -R~ l J . There is little difference between the profiles 
calculated at low and high redshift; this could in part be at- 
tributed to the cancellation of the nearly opposite effects of 
increasing halo concentration and decreasing mean halo den- 
sity with decreasing redshift. Galactic stellar mass loss from 
stellar evolution leads to a decrease of stellar surface density 
at all radii but does not modify the overall character of the 
final surface density profile. 

In Figure [2] left panel, we show the evolution of the NSC 
mass proxy Mio pc , here simply calculated as the stellar mass 
contained in a sphere with radius r = lOpc as a function of 
the initial migration time of t he A:th cluster to migrate to the 
galaxy center (see Section 3.3 I. The latter is an approximation 
since the migration time evolves as the cluster loses mass and 
as its orbit decays, and as the galactic potential is modified by 
the previous k— 1 migrating clusters. The initial value of Mio pc 



in Figure|2] ranging from ~ 10 2 M Q to > 1O 5 M , corresponds 
to the initial mass of the fastest-migrating cluster in each sim- 
ulation; forM max = 10 7 M Q , the mass of the fastest-migrating 
cluster to reach r < lOpc already falls within the range of ob- 
served NSC masses (see Section 5.1 1. The mass M^pc con- 
tains 96% of the mass of the unresolved nucleus M cent ac- 
counting for the cluster stars that have reached r = 1 pc and 
an additional "resolved" mass deposited by clusters migrating 
from lOpc to lpc. The mass Mio pc increases as an approxi- 
mate power-law in time with slopes between Mio pc oc t 1 ' 2 and 
M\o pc oc t and then levels off. In Table [T] we report the final 
asymptotic Mio pc and its ratio to the totalstellar mass, as well 
as the total time tin that it takes for Mio pc to reach one-half 
of its final value. In all the models presented in this work, 
the contribution of the unresolved central point mass M cent to 
M IOpc increases from about 40% for M max = 10 4 M Q to 86% 
for M max = 1O 7 M in spheroidals and from about 20% for 
M max = 10 4 M Q to 80% for M max = 1O 7 M in disks. We find 
that > 90% of the NSC mass is composed of clusters with ini- 
tial masses > 0.1M max , indicating that the lowest-mass clus- 
ters make up only a small fraction of the NSC mass and, thus, 
that the final central mass is not sensitive to M m ; n . 

The asymptotic M 1( ) pc is a strong function of the ICMF trun- 
cation mass scale M max , as is the time scale for approach to 
this asymptotic mass. Larger M max yield larger final central 
masses that are assembled on longer time scales; this is be- 
cause more massive clusters are able to migrate from larger 
radii. The central mass-to-ICMF truncation mass ratios are 
Mio pc /M max ~ 1— 1.5dex, and thus the central stellar system 
is assembled from many migrating clusters. The scaling of the 
resulting NSC mass with the ICMF characteristic mass scale 
is consistent with the trend recognized by Bekki (2010), who 
carried out -body integrations with an initial lognormal clus- 
ter mass function that approximates the evolved present-day 
globular cluster mass distribution. 

4.2. Disk Galaxies 

In the disk galaxy calculation we choose n = 1 and as- 
sume an approximate gas or stellar velocity dispersion in 
equation pi of a(r) = 0.2 [GM ga i(r)/r] 1/2 where M gal (r) = 
max[M+. $,M^ k (r)]+M DM (r) + M cea a(r) (see Section [33]). 
Figure [3] in its left and middle panels shows the projected 
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FIG. 2. — Total mass Mio pc accumulated after migration and disruption of star clusters within the central lOpc of a high-redshift spheroidal galaxy (left panel) 
and disk galaxy (right panel) as a function of migration time for those clusters that reach the central lOpc for M max = 10 4 Mq (solid line), M max = 10 s Mq (dashed 
line), 10 6 Mq (dash-dotted line), and 10 7 Mq (dotted line). The mass histo ries end when the last cl uste r makes it to the ce ntral lOpc. The shaded region in the 
left panel shows the lcr-range of NS C mass-to-total stel lar mass ratios from Ferrarese et al. (2006b I and Cote et al. (2006 1. In the right panel, the shaded region 
shows the range of NSC masses in Walcher et al.|j2005) . 




FIG. 3. — Same as Figure[T] but for our disk galaxy model at high redshift (left panel) and low redshift (middle panel), and at high redshift with 50% stellar 
mass loss (right panel). 



surface density profiles for stars (Mdi s k = 5 x 1O 9 M , R s = 
1.5kpc, n = 1), dark matter (M ha io = 5 x 10 10 M Q , high-z: c = 
3.5, r v j r = 20 kpc; low-z: c = 15, r v ; r = 100 kpc), and stars after 
migration and disruption for M max = (10 4 , 10 5 , 10 6 , 1O 7 )M . 
As we did for spheroidal galaxies in Section |4~T| here we test 
the effect of galactic mass loss due to stellar evolution in Fig- 
ure J3] right panel. The character of a disk galaxy's morpho- 
logical transformation in consequence of cluster migration is 
very similar to that observed in spheroidal galaxies. The den- 
sities at smaller radii < lOpc in disk galaxy models are some- 
what higher than in spheroidal galaxy models. This may be 
explained by a higher initial central stellar surface density in 
our disk models. In Figure|2| right panel, we show the evolu- 
tion of the mass enclosed within r = lOpc and in Table[2]we 
show the final, asymptotic mass Mio pc and the time scale t\/2 
for accumulation of a half the final mass in disk galaxy mod- 
els. We find that M 10pc /M max — 1 .7dex and that > 90% of the 



NSC mass is composed of clusters with initial masses above 
0.1M max forM raax > 1O 4 M indicating that, again, the central 
stellar system is assembled from many migrating clusters and 
that the NSC mass is not sensitive to M ra j n . 

4.3. Prompt Cluster Dissolution 

Observations of star cluster populations in nearby disk 
galaxies suggest that a large (~ 90%), though still highly un- 
certain, fraction of young clusters dissolve very quickly after 
formation, likely as the result of an initial loosely bound state 
and rapid expulsion of the natal gas cloud from the cluster 
by photoevaporation and radiation pressure (this is often re- 
ferred to as "infant mortality," see references cited in Section 
3.3 above). There is evidence that the likelihood of prompt 
dissolution is independent of cluster mass. To explore the im- 
pact of prompt dissolution on the morphological transforma- 
tion of the galaxy, we repeated the calculations carried out in 
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FIG. 4. — Same as Figure[T] but for a high-redshift proto-spheroidal galaxy (left panel) and proto-disk galaxy (right panel) with prompt cluster dissolution. In 
these models 90% of the clusters dissolve immediately independent of their mass, and the remaining 10% migrate inward prior to disruption. 
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FIG. 5. — The same as Figure|2] except that 90% of the clusters are dissolved immediately independent of their mass (see Section [4~3) . 



TABLE 1 

Spheroidal Galaxy" NSC Properties 



TABLE 2 
Disk Galaxy" NSC Properties 



Model 


log(M max ) 
(M ) 


l°g0l/2) 

(yr) 


log(M 10pc ) M 10pc /M sph 
(M ) (dex) 


MlOpc A/max 
(dex) 


Model 


log(M max ) 
(Mo) 


l°g(fl/2) 

(yr) 


log(M 10pc ) 


Wi 0pc /M disk 
(dex) 


M l0pc /M max 
(dex) 


High-z 


4 
5 
6 
7 


7.96 
8.60 
9.50 
9.34 


4.87 
6.32 
7.52 
8.20 


-4.13 
-2.68 
-1.48 
-0.80 


0.87 
1.32 
1.52 
1.20 


High-z 


4 
5 
6 
7 


7.49 
8.08 
8.87 
9.29 


5.67 
6.60 
7.77 
8.67 


-4.03 
-3.10 
-1.93 
-1.03 


1.67 
1.60 
1.77 
1.67 


Low-z 


4 
5 
6 
7 


8.09 
8.83 
9.59 
9.13 


5.24 
6.51 
7.64 
8.21 


-3.76 
-2.49 
-1.36 
-0.79 


1.24 
1.51 
1.64 
1.21 


Low-z 


4 
5 
6 
7 


7.44 
8.01 
8.92 
9.10 


5.75 
6.68 
7.86 
8.69 


-3.95 
-3.02 
-1.84 
-1.00 


1.75 
1.68 
1.86 
1.69 


a M sph = 


lO 9 M ,M ha i o 


= 10 10 M Q , 


no galactic mass 


loss or prompt dissolution. 


a M dlsk = 


5 X W 9 Mq, M h . Ao = 5 X 1O'°M , no 


galactic mass 


loss or prompt 



Sections |4.1| and |4.2| but assuming that 90% of randomly- 
selected clusters dissolve instantaneously and deposit their 
entire mass at the radius at which they formed. If prompt 
dissolution is mass-dependent with, for example, significantly 
less than 90% of high-mass clusters being instantaneously dis- 
solved, then we expect the resulting NSC mass from our sim- 



dissolution. 

ulations to be larger. The resulting stellar surface density pro- 
files are shown in Figure [4] Now, even for our largest ICMF 
truncation mass scale of M max = 10 7 M Q , visible evolution of 
the stellar surface density profile is confined within the inner 
~ 10-50pc. The time evolution of the mass contained in the 
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TABLE 3 

Spheroidal Galaxy" NSC Properties - 90% Prompt Dissolution 



TABLE 4 

Disk Galaxy 3 NSC Properties - 90% Prompt Dissolution 



Model 


log(M m ax) 
(M ) 


l°g('l/2) 

(yr) 


log(M 10pc ) 
(M & ) 


M 10p c/M sph 
(dex) 


Mio P c/M max 
(dex) 


Model 


log(M max ) 
(A/©) 


l°g0l/2) 

(yr) 


log(M 10pc ) 

(Me) 


(dex) 


sk A^opc/Mmax 
(dex) 


High-z 


4 
5 
6 
7 


6.66 
8.56 
9.94 
9.33 


4.26 
5.60 
6.55 
7.19 


-4.74 
-3.40 
-2.45 
-1.81 


0.26 
0.60 
0.55 
0.19 


High-z 


4 
5 
6 
7 


7.42 
8.20 
8.96 
9.75 


5.48 
5.92 
7.12 
7.75 


-4.22 
-3.78 
-2.58 
-1.95 


1.48 
0.92 
1.12 
0.75 


Low-z 


4 
5 
6 
7 


8.23 
9.14 
10.00 
9.63 


4.70 
5.87 
6.75 
7.24 


-4.30 
-3.13 
-2.25 
-1.76 


0.70 
0.87 
0.75 
0.24 


Low-z 


4 
5 
6 
7 


7.31 
8.25 
9.39 
9.36 


5.51 
6.05 
7.23 
7.85 


-4.19 
-3.65 
-2.47 
-1.85 


1.51 
1.05 
1.23 
0.85 


a M sph = 


10 9 Mq, Mhalo 


= 10 10 Mq, no galactic mass loss. 




" M disk = 


5x 10 9 M©,M halo = 5x 10'°M Q ,no 


galactic mass 


loss. 



innermost r = lOpc is shown in Figure[5] The central masses 
rise as Mio pc oc t 2 / 5 in most cases, and reach smaller asymp- 
totic values than without prompt dissolution. For M max = 
10 4 M Q there is no substantial departure from the inward- 
extrapolated Sersic profile. For M max = (10 5 - 10 7 )M o , the 
central masses decrease from ~ (10% -20%) of the values 
calculated without prompt dissolution in spherical and disk 
galaxies. In Tables [3jand [4|we summarize the properties of the 
NSCs (see SectionBT) in the simulations with 90% prompt 
cluster dissolution and find that Mio pc /M mllx ~ 0.2 -0.8 dex 
in spheroidals and Mio pc /M max ~ 0.8 -0.5 dex in disks. The 
timescale t t /2 for half of the final NSC mass to accumulate 
extends up to lOGyr. 

5. DISCUSSION 

5.1. Comparison with Observed NSCs and Pseudobulges 

Our model calculations show that cluster migration can 
bring about galactic morphological transformation in which 
a new stellar density component grows at the center of 
the galaxy. The new component is a drastic departure 
from the inward-extrapolated outer surface density profile 
of the galaxy, and this calls for comparison with NSCs in 
spheroidals and late-type disks, and with central light ex- 
cesses in spheroidals and pseudobulges in disks; we refer to 
both of the latter phenomena as "pseudobulges." The surface 
density profiles shown in Figures [T] [3] and|4]do not suggest 
a clear separation into a compact component that would be 
compared with an NSC and a more extended component that 
would be compared to a pseudobulge. We caution against 
overinterpreting the detailed profile because our model for 
progressive cluster mass loss and dissolution is crude and in- 
evitably fails to accurately account for the evolution of the 
internal structure of the cluster. For example, if the cluster 
is dense enough for its core to collapse through mass segre- 
gation and two-body relaxation, then the cluster core could 
maintain integrity longer and would migrate to smaller radii 
than in our calculations. More accurate estimates of the de- 
tailed innermost stellar density profile can only be achieved 
with A^-body simulations resolving the internal structure of 



the mig rating clusters (see, e.g., Capuzzo-Dolcetta & Miocchi 
2008a b). With this caveat in mind, we pursue comparison of 



only the more robust characteristics of the calculated profiles 
with those of observed galaxies. 

In Figures [2] and H] which show the evolution of the mass 
Mjopc contained in the innermost ten parsecs, we indicate the 
range of NSC masses corresponding to the fractions of the 
spheroidal luminosit y contained in the NS Cs i n the Virgo 
Cluste r Survey of |Ferrarese et al. ( 2006b I and |C6te et aL] 
( |2006[ ), and also indicate the range oFabsolute NSC masses 



in the survey o f late -type disks by Bo ker et al.| (12004} and 
( |2005[l . If Miopc is indeed a valid proxy 



Walcher et al.l 



tor NSC mass, then this allows us to identify the range of 
ICMF truncation mass scales M max consistent with the ob- 
served NSCs. The models not allowing for prompt dissolu- 
tion are consistent with M max ~ 10 4 - 1O 5 M for spheroidals 
and M max - 10 4 - 10 6 M Q for disks. The models with 90% 
prompt dissolution are consistent with M max ~ 10 5 - 10 6 M Q 
for spheroidals and M max ~ 10 5 - 10 7 M Q for disks. These es- 
timates are consistent with the observationally inferred and 
theoretically a ntici pated values of the ICMF truncation mass 
scale (Section |2.1[ i. The results uniformly exclude the possi- 
bility that the ICMF truncation mass scale in spheroidals and 
late-type disks is above 1O 7 M0. 

The surface density profile modified by cluster migration 
already starts departing upward from the inward-extrapolated 
outer Sersic (or exponential) law at a radius that increases 
with M max = (10 4 - 10 6 )M Q from - 20pc to - 300pc in the 
models without galactic mass loss and prompt dissolution. 
With galactic mass loss, the departure radius ranges from 
< lOpc to ~ 200pc for the same range of M max . These depar- 
ture radii are relatively small fractions of the disk exponen- 
tial scale length R s in disk galaxies. In the models with 90% 
prompt dissolution, the departure radius ranges from ~ lOpc 
to - 40pc for M max = (10 5 - 1O 7 )M . For M max = 1O 4 M , no 
significant surface density excess is present. Because of the 
crudeness of the prescriptions that we employ to model cluster 
migration and dissolution, and the sensitivity of the innermost 
density profile to these prescriptions, we do not attempt to fit 
an analytic profile, such as a two-component Sersic profile al- 
lowing for a photometrically distinct central light excess (and 
perhaps a third component-the NSC), to the surface density 
profile of the final galaxy. Nevertheless, the excess we ob- 
serve in the central one or few hundred parsecs is suggestive 
of pseudobulge s that have been identified photometrically in 
disk galaxies (Kormendy & Kennicutt 2004; Fisher & Drory 
12008] |Fisher et al.|2009[|Fisher & Drory|2010||Weinzirl et al. 
2009 ) and the "central light excesses" identified in bulgeless 
disks (Bok er et al.||2003[), incl uding the late-type-disk M33 
HKent|1987[|Minniti et al.|T993] >. 

The effective radii of the pseudobulge components identi- 
fied in the recent surveys by Fisher et al. and Weinzirl et al. 
seem to be compatible with our more optimistic models that 
ignore prompt dissolution, and are on average larger than the 
radii within which we detect surface density excess in the pes- 
simistic models with 90% prompt dissolution. We caution 
against direct comparison because in the present work, in an 
attempt to emphasize sensitivity to the variation of the ICMF 
truncation mass scale M max , we have held the parameters of 
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our dark halo and initial baryonic disk (or spheroid) fixed at 
values that seem to correspond to galaxies that are somewhat 
smaller than the typical pseudobulge hosts. We can only con- 
clude that a pseudobulge-like central stellar surface density 
increase is generic and that cluster migration is one potential 
contributor to pseudobulge assembly in disk galaxies, while 
other processes, such as angular momentum transport by stel- 
lar and gaseous bars, certainly also contribute, in line with 
the observation that pseudobulge hosts g enerally have nuclear 
bars, rings, or nuclear spi rals (e.g., |Kormendy & Kennicutt| 
2004, Fis her & Drory|2008) . 

The most massive clusters that we have considered are 
still substantially less massive than the giant ~ 10 8 - 10 9 M Q 
clumps that are observed to be present and are theoreti- 
cally expected to be forming in globally gravitationally un- 
stable, rapidly-star-forming massive disks at high redshift 
(e.g., |Noguchi 1999; Bournaud et al. 2007; Elmegree n et al.| 
[2008i|Dekel et al.|2009b||Tacconi et al.|2010| and references 
therein). The super star clusters forming in these giant clumps 
should be more immune to dissolution in the tidal field of the 
galaxy and could reach the galactic central region intact. We 
speculate that there could be a critical characteristic ICMF 
mass scale above which clusters migrate intact and merge 
to produce a classical bulge (see, e.g., |Immeli et al.||2004t 



Elmegreen et al. 2008 Ceverino et al. 2010), and below which 
they suffer substantial mass loss en route to the galactic center 
and thus give rise to a pseudobulge. 



5.2. Implications for Massive Black Holes 
The apparent agreement of NSC-mas s-to-galactic-stellar 



mass ratios in spheroidals (~ 2 x 10 3 ; Cote et al.|200"6 



Fer- 



|rarese et al.|[20 06b; We hner & Harris||2006| l and massive- 
black-hole-to-galactic stellar mass ratios in ellipticals and 
bulges (~ [1-2] x 10" 3 ; e.g., |Kormendy & Richstone||l995l 



Wand el|[T999| Kormendy & Gebhardtp OOl ; Merritt & Fer- 
rarese|2001||McLure & Dunlop|2002[|Marconi & Hunt|2003[ ' 
Haring & Rix 2004 ) has prompted speculation that the same 



process may be responsible for the formation of NSCs and 
black holes. While the formation and growth of a massive 
black hole undoubtedly requires a gas-dynamical, dissipative 
process, our results suggest that an NSC can be assembled 
nondissipatively, by the collisionless migration of star clus- 
ters, and thus, the observed agreement could be a coincidence. 

Although the well-studied NSC host galaxy M3 3 does not 
contain a central ma ssive black hole (Me rritt et al.||200"T) 
|Gebhardt et al.||2001| and references th erein), another one, 
NGC4395, doeslFiTi ppenko & Ho|2003[l , and still others con- 
tain AGNs ( |Seth et al.|2008| ). Since AGNs and the growth of 
a central black hole require gas inflow into the center of the 
galaxy, NSC growth from migrating disk clusters would not 
generally be accompanied with black hole growth, because 
disk clusters contribute stellar mass without augmenting the 
black hole mass (although gas in flow may be enhanced by 
the migrating clusters, see, e.g., Goodman & Rafikov||200l( 
Chang 2008). This suggests that the central black hole mass 
in bulgeless disks should not be correlated with the mass 
of the NSC and that any pseudobulge that is present, un- 
less the migrating clusters independently synthesize massive 
(or intermediate-mass) black holes which, in this case, they 
would deliver to t he center to merge to fo rm a more massive 
central black hole (Elmeg reen et al.|2008| l. 



6. CONCLUSIONS 

We studied the morphological transformation of isolated 
spheroidal and late-type disk galaxies driven by migration and 
tidal disruption of embedded star clusters. With the help of 
empirically and theoretically calibrated cluster migration and 
disruption time scales, we have tracked the mass that accu- 
mulates in the central region of a spheroidal or late-type disk 
galaxy and have compared the resulting innermost surface 
density profile to that of NSC and pseudobulge host galax- 
ies. We have focused on the variation in the degree of galaxy 
transformation with the ICMF truncation mass scale; the as- 
sumptions and limitations of our model are described in Sec- 
tion !3.4l Our main conclusions are: 

1. The amount stellar mass that the migrating clusters 
transport into the central few tens or hundreds of parsecs is 
sensitive to the maximum mass of the clusters forming in 
the galaxy. Larger ICMF high-mass truncation mass, M max , 
scales yield larger and more spatially extended central stellar 
accumulations. 

2. Because our model for progressive cluster disruption is 
crude, we cannot attach significance to the detailed shape of 
the innermost surface density profile (the number of photo- 
metrically distinct components, inflection points, etc.), but we 
do compare the robust features of our synthetic profiles with 
the observed properties of NSC host galaxies. 

3. Our model yields NSC masses compatible with the ob- 
served NSC masses and scaling relations if the ICMF trunca- 
tion mass scale is 10 5 M Q < M max < 1O 6 M , which is con- 
sistent with direct M max estimates in the nearby star-forming 
galaxies. To match NSC masses in galaxies with a higher de- 
gree of prompt dissolution ("infant mortality"), higher ICMF 
truncation mass scales are needed. 

4. The surface density profiles modified by cluster migra- 
tion exhibit excesses above the inward-extrapolated outer Ser- 
sic (or exponential) profile in the central few tens to few hun- 
dreds of parsecs. These excesses are suggestive of the pseu- 
dobulge phenomenon in disk galaxies. The formation of NSCs 
and central surface density excesses in bulgeless galaxies is 
inevitable and generic simply as a consequence of the clus- 
tered nature of star formation. 

We thank D. Fisher, M. Gieles, J. Kormendy, and C. J. 
Walcher for invaluable discussions and also thank the organiz- 
ers of the workshop Nuclear Star Clusters across the Hubble 
Sequence at MPIA, Heidelberg, 2008. We would also like to 
thank an anonymous referee for helpful comments on an ear- 
lier version of this work. We acknowledge NSF grant AST- 
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